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Epigenetic dysfunction is a known contributor in carcinogenesis, and is emerging as a mechanism involved in toxicant- 
induced malignant transformation for environmental carcinogens such as arsenicals or cadmium. In addition to aberrant 
DNA methylation of single genes, another manifestation of epigenetic dysfunction in cancer is agglomerative DNA 
methylation, which can participate in long-range epigenetic silencing that targets many neighboring genes and has 
been shown to occur in several types of clinical cancers. Using in vitro model systems of toxicant-induced malignant 
transformation, we found hundreds of aberrant DNA methylation events that emerge during malignant transformation, 
some of which occur in an agglomerative fashion. In an arsenite-transformed prostate epithelial cell line, the protocadherin 
(PCDH), HOXC and HOXD gene family clusters are targeted for agglomerative DNA methylation. The agglomerative 
DNA methylation changes induced by arsenicals appear to be common and clinically relevant events, since they occur 
in other human cancer cell lines and models of malignant transformation, as well as clinical cancer specimens. Aberrant 
DNA methylation in general occurred more often within histone H3 lysine-27 trimethylation stem cell domains. We found 
a striking association between enrichment of histone H3 lysine-9 trimethylation stem cell domains and toxicant-induced 
agglomerative DNA methylation, suggesting these epigenetic modifications may become aberrantly linked during 
malignant transformation. In summary, we found an association between toxicant-induced malignant transformation 
and agglomerative DNA methylation, which lends further support to the hypothesis that epigenetic dysfunction plays an 
important role in toxicant-induced malignant transformation. 



Introduction 

Arsenic and cadmium are known carcinogens that are associated 
with many other chronic diseases; unsafe exposure levels to these 
toxicants remain a problem for millions of people throughout 
the world. 1 " 4 Experiments have shown that long-term arsenic or 
cadmium exposure in vitro induces malignant transformation of 
several immortal cell lines, but the underlying mechanisms driv- 
ing these toxicant-induced malignant transformations are multi- 
faceted and remain elusive. 5 " 8 Epigenetic dysfunction is emerging 
as an important aspect of both arsenic and cadmium associated 
chronic diseases, especially cancer. 9 " 11 

In addition to genetic alterations, epigenetic changes are a 
driving force in cancer. 12 Earlier studies have shown that the aber- 
rant DNA methylation profiles of distinct cancer models correlate 
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with the degrees of their aberrant phenotypes. 13,14 In other words, 
as the phenotype of a cell progresses toward malignancy, its epi- 
genetic patterning also becomes increasingly more distinct from 
that of its normal tissue counterpart. For example, during the 
immortalization of human mammary epithelial cells few to 
hundreds of differentially methylated regions (DMRs) arose 
concurrently and were correlated with the escape from distinct 
proliferation barriers. 14 During the later stages of tumorigenic 
progression, DMRs were observed to accumulate more gradually 
while correlating with advanced malignant phenotypes. 13,14 

In addition to targeting individual genes, cancer-associated 
epigenetic dysfunction is known to occur across large chromo- 
somal regions which can encompass many genes. In normal 
cells, epigenetic mechanisms such as DNA methylation, his- 
tone H3 lysine-27 trimethylation (H3K27me3) and histone 
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Table 1. Cell lines and samples analyzed by MeDIP-on-Chip 



Sample 


Exposure 


Cone. 


Duration 


Culture post-exposure 


Tumorigenic 1 


Reference 


HUC (n = 2) 


None 


NA 


NA 


NA 


ND 




UROtsa (n = 3) 


None 


NA 


NA 


NA 


No 


25, 55 


URO-MSC12+24(-) 


MMA 1 " 


50 nM 


12 weeks 


24 weeks 


Yes 


31 


URO-MSC52 


MMA 1 " 


50 nM 


52 weeks 


Yes 


Yes 


56 


URO-ASSC 


As 111 


1 u.M 


52 weeks 


Yes 


Yes 


5 


URO-CDSC 


Cd" 


1 |j,M 


52 weeks 


Yes 


Yes 


5 


BFTC-905 


As 111 


NA 


NA 


Yes 


Yes 


39 


BFTC-909 


As 111 


NA 


NA 


Yes 


Yes 


39 


Bladder tumor (n = 6) 


NA 


NA 


NA 


NA 


Yes 




PrEC (n = 2) 


None 


NA 


NA 


NA 


ND 




RWPE-1 (n = 2) 


None 


NA 


NA 


NA 


No 


26 


CAsE-PE (n = 2) 


As 111 


5 |xM 


33 weeks 


None 


Yes 


6,35 


CTPE 


Cd" 


10 (jtM 


10 weeks 


None 


Yes 


8,57 


WPE1-NB26 


MNU 


100 (jug/mL 


4 cycles of 1 h 


Yes 


Yes 


51 


LNCaP 


None 


NA 


NA 


NA 


Yes 


58 


PC3 


None 


NA 


NA 


NA 


Yes 


59 


DU145 


None 


NA 


NA 


NA 


Yes 


60 



NA, not applicable; ND, no data. a Tumorigenic determination for the UROtsa and RWPE-1 cell lines including their derived counterparts is based on 
their ability to form tumors in immunocompromised mice. The remaining tumorigenic samples are listed as such because each was derived or isolated 
from a human tumor. 



H3 lysine-9 trimethylation (H3K9me3) play important roles in 
silencing individual genes as well as large chromosomal regions 
or even entire chromosomes such as the inactive X of females 
(reviewed in ref. 15). Likewise, genomic scale analyses of cancers 
have shown that epigenetic dysfunction targets individual genes 
and larger chromosomal regions containing many genes. 16 " 20 
These large scale aberrations of repressive epigenetic modifica- 
tions have collectively been called long-range epigenetic silenc- 
ing (LRES). 21 Groups of DMRs that are enriched within defined 
chromosomal regions have been called agglomerative epigenetic 
aberrations and these agglomerative DNA methylation events 
can contribute to LRES. 19 LRES events are unique epigenetic 
lesions that have been observed in various types of tumors, 
but they have not been linked to any specific carcinogen expo- 
sure 16,19,22-24 »p ne £ act jhgy are ODSer ved in tumors supports 

their clinical relevancy and suggests that carcinogen exposure 
could lead to their induction. 

Our work examines two in vitro models, RWPE-1 and 
UROtsa, which have been used successfully to study epigen- 
etic dysfunction associated with arsenical or cadmium-induced 
malignant transformation. These immortalized cell lines are 
non-tumorigenic and were derived from the prostate and ure- 
ter respectively, which are target sites for cadmium and arsenic- 
induced cancers. 25 " 30 Long-term, independent exposures of these 
cell lines to arsenite (As 1 "), monomethylarsonous acid (MMA 111 ) 
and cadmium (Cd 11 ) induced hyperproliferation, conferred 
anchorage independent growth and tumorigenicity. 5 ' 6 ' 8,31,32 The 
resulting tumorigenic cell lines provide valuable models to study 
the epigenetic component of arsenical or cadmium-induced 
malignant transformation in relevant target tissues. Previous 



studies have shown that the arsenical or cadmium-transformed 
variants of RWPE-1 and UROtsa harbor stable aberrations to 
their DNA methylation profiles, including global hypomethyl- 
ation and gene promoter hypermethylation. 13 ' 33 " 35 

We extend prior studies using higher resolution DNA methyl- 
ation analyses, allowing us to gain further insights into aberrant 
DNA methylation patterns induced by long-term toxicant expo- 
sure. In this study we compare aberrant DNA methylation in 
two relevant models of arsenical or cadmium-induced malignant 
transformation, which enables us to distinguish the effects which 
are model-specific, toxicant-specific or generally found in can- 
cers. Our results lend further support to the hypothesis that epi- 
genetic dysfunction plays an important role in toxicant-induced 
malignant transformation by showing that clinically relevant, 
agglomerative aberrations of DNA methylation occur during 
arsenical and cadmium-induced malignant transformation. 

Results 

Malignant transformation by toxicants induces hundreds of 
DMRs in distinct models. This study was an effort to better 
understand the DNA methylation events of arsenic and cad- 
mium-induced malignant transformation. Using MeDIP-on- 
Chip we analyzed the DNA methylation changes that occur 
during arsenic or cadmium-induced malignant transformation 
of immortal cell lines. Table 1 lists the cell lines that were ana- 
lyzed by MeDIP-on-Chip in this study with descriptions of their 
treatments and phenotypes. The schematic in Figure 1A shows 
the lineage of the in vitro models of toxicant target tissues used 
in this study which represent normal, immortal and malignant 
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phenotypes. HUC and PrEC are pheno- 
typically normal cells with finite lifes- 
pans. Although these cell strains have 
distinct genotypes, their DNA meth- 
ylation profiles can be used to represent 
that of the normal epithelial tissues of 
origin of UROtsa and RWPE-1 due to 
the stability and low variability of DNA 
methylation in normal tissues. 36,37 We 
used this collection of normal, immor- 
tal and malignant models to help us put 
toxicant-induced aberrant DNA meth- 
ylation into context with the meth- 
ylation changes that occur generally in 
prostate and urothelial carcinogenesis. 

To understand the context of the 
collective DNA methylation change 
associated with arsenical or cadmium- 
induced malignant transformation, we 
acquired a genome-wide DNA methyla- 
tion profile for each of the cell lines and 
bladder tumors using MeDIPs coupled 
to 470,000 probe tiling human pro- 
moter microarrays that query 17,000 
human gene promoters. We then used 
multidimensional scaling (MDS) to 
summarize the data from individual 
microarray probes and thus graphi- 
cally represent the variability in the 
DNA methylation profiles between 
all the samples (Fig. IB). The finite- 
lifespan cells, HUC and PrEC, are in 
close proximity to one another after 
MDS, which indicates that they have 
similar DNA methylation profiles. 
MDS also shows that the immortal- 
ized RWPE-1 and UROtsa cell lines do 
not cluster with the finite-lifespan cells, 
indicating that immortalization had a 
pronounced effect on the DNA meth- 
ylation profiles of these cell lines, which 
is consistent with earlier studies that 
describe aberrant DNA methylation 
associated with immortalization. 14,38 Of 
the two immortalized models, RWPE-1 
has a methylation profile more simi- 
lar to that of the finite-lifespan cells. 
In summary, the MDS shows that the 
normal cells have similar methylation 
profiles while the distinct immortal- 
ized cell lines have variable degrees of 
aberrant DNA methylation, which 
may be related to their pathways of 
immortalization. 

Our analysis includes five cancer 
cell lines that collectively represent the 
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Figure 1. (A) Schematic showing the origins of the toxicant-transformed cell lines. The dashed 
arrows indicate that the finite-lifespan and immortal cells have distinct genotypes whereas the 
malignant cell lines are direct descendants of their corresponding immortal parent. The immortal 
RWPE-1 cell line was created from a primary culture of PrEC that was immortalized by the addition 
of HPV18. RWPE-1 cultures were chronically exposed to arsenite or cadmium to create the malignant 
cell lines CAsE-PE and CTPE. WPE1-NB26 cell were created by 100 (jig/mL MNU treatment. The im- 
mortal UROtsa cell line was created from a primary culture of HUC by the addition of an SV40 large-T 
antigen construct. UROtsa cells were malignantly transformed via long-term exposures to arsenite, 
MMA 1 " and cadmium giving rise to URO-ASSC, URO-MSC and URO-CDSC respectively. (B) Multidi- 
mensional scaling of pairwise distances derived from the DNA methylation level of gene promoters 
of all samples shows the relative similarity of the overall DNA methylation profiles of all the samples. 
The distances between samples in the MDS plot represent the degree of difference between the 
samples' methylation profiles. The six bladder tumor biopsies are labeled as CD450, CD409, CD293, 
CD196, CD095 and CD428. Arsenical or cadmium-induced malignant transformation of RWPE-1 
alters the DNA methylation profile so that it more closely resembles the cancer cell lines and tumor 
biopsies. 
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aberrant DNA methylation events that occur in carcinogenesis. 
PC3, LNCaP and DU145 were each derived from metastatic 
prostate cancers and as such their methylation profiles reflect 
metastatic prostate cancer. The BFTC cell lines were derived 
from urothelial carcinomas of individuals who lived in Blackfoot 
disease endemic regions, which are regions specifically linked to 
high levels of environmental arsenic. 39 The prostate cancer cell 
lines and BFTC cell lines have the same tissue origins as RWPE-1 
and UROtsa, and thus represent tissue matched models of the 
end stages of prostate and urothelial carcinogenesis. These five 
cancer cell lines are closer to one another than to the immortal 
or finite-lifespan cells in the MDS plot, which suggests that a 
portion of their DNA methylation profiles is common between 
them, but different from the other stages (e.g., finite-lifespan, 
immortal) (Fig. IB). This is confirmed in Figure 2 which shows 
that the BFTC cell lines and the prostate cancer cell lines have a 
substantial overlap of differentially methylated probes compared 
with their respective finite-lifespan cells and 30% of these probes 
are not differentially methylated in immortalized cells compared 
with normal cells (Fig. SI). 

In order to assure that the differences between the finite- 
lifespan cells and the cancer cell lines are relevant to clinical can- 
cers, we included six human bladder tumor biopsies in our analy- 
sis. After MDS, the six tumors occupy a common region of the 
plot (Fig. IB). This distribution of the tumor biopsies indicates 
that their DNA methylation profiles vary from one another while 
maintaining some similarity relative to the other sample types. 
Overall, the tumor biopsies are nearer to the cancer cell lines 
than to the finite lifespan cells which suggests that there may 
be a common set of aberrant DNA methylation events shared by 
both the clinical samples and cancer cell line models that corre- 
lates with their malignant phenotypes. Venn diagrams compar- 
ing aberrantly methylated probes between groups show that 26% 
of the differentially methylated probes in the tumor biopsy group 
are also differentially methylated in BFTC and prostate cancer 
cell lines, and 98% of these probes are hypermethylated (Fig. 2). 

A comparison of the DNA methylation profiles of the virally- 
immortalized, non-tumorigenic cell lines RWPE-1 and UROtsa 
with their toxicant-transformed counterparts affords an estimate 
of the degree of aberrant DNA methylation induced by each 
toxicant. The MDS shows that in every case, transformation 
by toxicant exposure affected the DNA methylation profile in a 
similar way for each cell type. Each of the toxicant-transformed 
cell lines shifts in a similar direction relative to their untreated 
counterparts (Fig. IB), which suggests that each has an analo- 
gous methylation change in a comparable portion of the probes. 
The DNA methylation changes that occur with toxicant-induced 
malignant transformation of RWPE-1 cells result in DNA meth- 
ylation profiles that are more like cancer cell lines and tumor 
biopsies than before toxicant-induced transformation (Fig. IB). 
The increasing similarity between DNA methylation profiles 
of toxicant-transformed RWPE-1 cells and those of cancer cell 
lines and tumor biopsies is driven mostly by hypermethylation 
since there are more hypermethylated probes and they overlap 
more often than hypomethylated probes (Fig. 2). Taken together, 
the MDS and Venn diagrams show that toxicants may induce 



Figure 2. Venn Diagrams showing the numbers of significant probes 
compared with respective parent or normal that are shared between 
the different groups of malignant samples. (A) Hypermethylated 
probes commonly overlap between groups, which indicates that there 
are common targets of hypermethylation in the malignant samples. 
(B) There are fewer hypomethylated probes overall and they overlap 
less frequently which indicates that hypomethylation has fewer com- 
mon targets between the malignant sample types. 
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aberrant DNA methylation events that are common in cancer. 
In summary, an important interval of epigenetic change can be 
attributed to toxicant-induced transformation, and this interval 
correlates with the acquisition of a malignant phenotype. 

To more finely resolve the epigenetic changes induced by toxi- 
cants and those associated with malignant transformation we 
identified specific DMRs between grouped samples. Compared 
with their parent cells, we found a total of 503 and 409 hyper 
DMRs in the CAsE-PE cells and arsenical-transformed UROtsa 
cells respectively (Fig. 3). Hypo DMRs were fewer, numbering 
72 in CAsE-PE and 98 in the arsenical-transformed UROtsa 
lines. We also analyzed one cadmium-transformed variant of 
each immortal cell line as well as the WPE1-NB26 cells which 
showed that the number of cadmium or MNU-induced DMRs 
is comparable to those induced by arsenic. In summary, toxicant- 
induced malignant transformation causes hundreds of DMRs in 
distinct target tissues. 

Agglomerates of aberrant DNA methylation in toxicant- 
transformed models. We next set out to locate extended regions 
of the genome that are enriched for DNA hypermethylation 
events after arsenical transformation of the RWPE-1 cells. This 
model was chosen to lead the study because the RWPE-1 cells 
had a relatively normal DNA methylation profile before expo- 
sure. In order to find agglomerative regions of aberrant DNA 
methylation we systematically scanned the MeDIP-on-Chip data 
by searching for genomic windows that contained more DMRs 
than expected by chance. Figure 4 shows in red the genomic 
locations where there is a significant enrichment of DMRs from 
the CAsE-PE cells. We found six genomic regions with signifi- 
cant enrichment of hyper DMRs in the CAsE-PE cells. These six 
genomic regions were further cropped to contain only runs of four 
or more DMRs with less than 200 kb between adjacent DMRs 
which left only three aberrant methylation events that meet these 
criteria (Table 2). These three agglomerative events target groups 
of consecutive genes in the protocadherin, HOXC and HOXD 
gene family clusters. The methylation status determined by 
MeDIP-on-Chip was confirmed for several of the DMRs within 
agglomerative regions using the MassARRAY EpiTYPER assay 
(Fig. 5; Figs. S2 and S3). The presence of these agglomerative 
events in the CAsE-PE cells indicates that arsenic-induced aber- 
rant DNA methylation can target large chromosomal regions in 
addition to individual genes during malignant transformation. 

Each of the three genomic regions with agglomerative hyper 
DMRs in CAsE-PE cells also shows aberrant hypermethylation 
in the BFTC cell lines, bladder tumor samples and prostate can- 
cer cell lines relative to HUC or PrEC (Table 2). These regions 
are also hypermethylated after cadmium or MNU-induced 
transformations of RWPE-1 (Table 2, CTPE and WPE1-NB26). 
The commonality of hypermethylation in these regions suggests 
that they are not targets unique to arsenic-transformed cells, but 
are hypermethylated generally in prostate and urothelial cancers. 
The fact that some of the agglomerative DMRs are also observed 
in the bladder tumor biopsies, prostate cancer cell lines, black- 
foot cell lines, CTPE and WPE1-NB26 cells suggests that they 
are clinically relevant, bona fide epigenetic lesions that can be 
induced by exposure to multiple distinct toxicants (Figs. S4-8). 



Figure 3. DMRs were calculated between the groups of samples shown. 
The blue and yellow bars show the number of hyper and hypo DMRs 
generated by the corresponding contrast. 



The greatest enrichment of DNA hypermethylation in 
CAsE-PE cells occurred at the PCDH gene family cluster on 
chromosome 5 and at the HOXC gene family cluster on chromo- 
some 12 (Fig. 4). The HOXC and PCDH gene family clusters 
are decisive examples of agglomerative, aberrant DNA methyla- 
tion. The PCDH cluster has previously been identified in breast 
cancer as a target of agglomerative DNA hypermethylation and 
in Wilms' and colorectal tumors as a target of long-range epi- 
genetic silencing. 19,23,24 The HOX clusters are also known targets 
of LRES, 19,20 and specifically HOXC hypermethylation occurs in 

in 

astrocytomas. 

To better characterize the degree of aberrant methylation 
within the PCDH, HOXC and HOXD gene family clusters 
we hybridized independent MeDIPs of our samples to an addi- 
tional Agilent tiling microarray of custom design. This array has 
extensive tiled coverage of all the HOX and PCDH gene family 
clusters giving more complete data with respect to the changing 
methylation pattern in these regions. On this custom array we 
found that in the CAsE-PE cells, the protocadherin G subfam- 
ily cluster (PCDHG) harbors 18 DMRs that collectively span 
16 kb (Fig. 5). In the HOXC region there are 23 DMRs which 
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Figure 4. Hyper DMRs from CAsE-PE cells were analyzed for agglomerates. The nega- 
tive log 10 of the p value calculated for each window was plotted as a distance from the 
respective chromosomal locations. The same calculation was applied to H3K9me3 
stem cell domains to find regions of enrichment within the Agilent promoter ChlP- 
on-Chip microarray coverage. Points were plotted in red above or green below the 
corresponding chromosomal location when the enrichment of DMRs or H3K9me3 
domains within the window was statistically significant (adjusted p < 0.05). 



collectively span 28.4 kb of the 151 kb total array coverage in 
this region (Fig. S2). Within these regions, there is at least one 
DMR for nearly every gene promoter and these DMRs associate 



with the CpG islands of the promoters. The HOXD 
region contains only five DMRs (Fig. S3). 

Association of arsenical-induced agglomera- 
tive DMRs with histone modification domains 
of embryonic stem cells. In order to better under- 
stand the potential mechanisms of arsenic-induced 
epigenetic dysfunction, we compared the DMRs 
found in CAsE-PE cells with known domains of his- 
tone modifications in human embryonic stem cells. 
H3K27me3 and H3K9me3 domains of stem cells 
were determined previously. 41 H3K27me3 is a tran- 
scriptionally repressive mark deposited by the poly- 
comb group proteins and genes that are polycomb 
targets in stem cells are often targets of aberrant 
DNA methylation in tumors. 42,43 H3K9me3 is also 
a repressive mark and an important component of 
heterochromatin which has been shown to work in 
concert with DNA methylation in part by directing 
its spreading. 44 " 46 Comparing the locations of DMRs 
to H3K27me3 and H3K9me3 stem cell domains 
may help determine if aberrant DNA methylation is 
linked to either of these histone modifications. 

To assess the spatial relationships between the 
arsenical-induced hypermethylation and the stem 
cell histone modification domains we calculated 
the number of hypermethylated microarray probes 
from CAsE-PE cells that occur within H3K27me3 
or H3K9me3 stem cell domains. Figure 6 shows that 
76 and 18 percent of the hypermethylated probes in 
CAsE-PE cells overlap H3K27me3 and H3K9me3 
stem cell domains, respectively. These overlaps are 
highly significant based on the expected overlaps of 
27 and four percent for H3K27me3 and H3K9me3 
domains, respectively. These results suggest that 
H3K27me3 and H3K9me3 stem cell domains are 
susceptible to aberrant DNA methylation during 
malignant transformation of CAsE-PE cells. 

Next we assessed whether the CAsE-PE agglom- 
erative DNA methylation events associate with 
genomic regions enriched with H3K27me3 or 
H3K9me3 stem cell domains. To do this we scanned 
the portion of the genome covered by the human 
promoter microarray for extended regions that 
were enriched with these domains and compared 
them to the locations of the agglomerative DMRs. 
Enriched regions of H3K27me3 domains partially 
overlapped the CAsE-PE agglomerative events in 
an expected manner considering that H3K27me3 
stem cell domains cover almost 27 percent of the 
microarray. There was a striking overlap of the 
H3K9me3 enriched regions with the agglomerative 
DNA methylation events. Figure 4 shows in green 
the genomic locations where there is a significant 
enrichment of H3K9me3 stem cell domains. In the CAsE-PE 
cells, all three of the agglomerative DNA methylation events 
overlap with significant peaks of H3K9me3 domain enrichment 
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Table 2. Genomic regions of arsenic-induced agglomerative DNA hypermethylation in CAsE-PE cells (n = 2) found with the Agilent promoter ChlP-on- 
Chip microarray 



Chr From (bp) 


To (bp) 


Size 
(kbp) 


#of 
DMRs 


Genes with DMRs 


Other groups or samples 
with hypermethylation 


chr2 176678716 


176760984 


82 


5 


HOXD11, HOXD10, HOXD3, HOXD1 


(LNCaP, PC3 and DU145), CTPE, (bladder 
tumors), [URO-ASSC, URO-MSC12+24H and 
URO-MSC52], URO-CDSC, (BFTC-905 and 



chr 5 139907571 140844320 



EIF4EBP3, CD14, PCDHA3, PCDHAC2, 
PCDHB6, PCDHB7, PCDHB16, PCDHGA2, 
937 18 PCDHGA3, PCDHGB4, PCDHGA8, 

PCDHGB5, PCDHGA10, PCDHGB7, 
PCDHGA11, PCDHGA1 2, PCDHGC4 



(LNCaP, PC3 and DU145), CTPE, (bladder 
tumors), (BFTC-905 and BFTC-909), 
WPE1-NB26 



chr 12 52616802 



HOXC13, HOXC12, HOXC11, HOXC10, 
52732432 116 15 hsa-mir-196a2, HOXC9, HOXC8, HOXC6, 

HOXC5, HOXC4 



(LNCaP, PC 3 and DU145), CTPE, (bladder 
tumors), [URO-ASSC, URO-MSC12+24H and 
URO-MSC52], URO-CDSC, (BFTC-905 and 
BFTC-909), WPE1-NB26 



Consecutive genomic windows with p < 0.05 (Fig. 4) were trimmed to include only the regions that contain runs of 4 or more DMRs with less than 
200 kb between adjacent DMRs. This filter was based on randomized CAsE-PE hyper DMR data, which showed that the expected number of runs with 
four or more DMRs and 200 kb or less between them is zero. In the non-randomized CAsE-PE data there are three runs that meet these criteria and oc- 
cur within significant windows of hyper DMR enrichment. Other samples with hyper DMRs in these regions are shown in the final column. Samples in 
parenthesis were analyzed as groups. 



(Fig. 4). These occur at the HOXD cluster on chromosome 2, 
the PCDH cluster on chromosome 5 and the HOXC cluster on 
chromosome 12. These overlaps suggest a potential connection 
between H3K9me3 stem cell domains and agglomerative DNA 
methylation. 

Discussion 

This study reports agglomerative DNA hypermethylation in 
toxicant-transformed cell lines, linking chronic arsenic and cad- 
mium exposures in vitro with hypermethylation of large chro- 
mosomal regions. This new association supports the hypothesis 
that epigenetic dysfunction plays an important role in arsenic- 
induced malignant transformation because the agglomerative 
events are clinically relevant, occurring in human tumor samples 
(Fig. S5). Agglomerative DNA hypermethylation is one modi- 
fication that is sometimes involved in LRES. 18 " 20 LRES events 
have been observed in various human tumors and they reduce 
the transcriptional plasticity of groups of consecutive genes. 47 
While our study does not specifically address transcriptional 
silencing caused by these arsenic-induced agglomerative events, 
it does identify regions that are known sites of LRES in other 
models. 19,20 ' 23 ' 24 We suspect that agglomerative DNA methyla- 
tion often occurs at genes whose expression is already repressed 
by some other mechanism. Based on this study, we propose 
that while the agglomerative DNA methylation events are not 
always driving gene expression changes, they may represent deep 
epigenomic scars left by malignancy driving catastrophic events 
caused by toxicant exposure. 

The agglomerative DNA hypermethylation events that occur 
at the PCDH and HOXC clusters were independently induced 
during arsenite, cadmium and MNU mediated malignant 



transformations, suggesting that they are reproducible events 
(Fig. 4; Figs. S4 and S8). Therefore, the RWPE-l/CAsE-PE 
cell line model can be used as a tractable system in which to 
study the mechanisms involved in agglomerative DNA methyla- 
tion and LRES, especially in the candidate regions found here 
such as the HOX and PCDH gene family clusters. The repro- 
ducibility of these events by distinct toxicants also suggests that 
each perturbs some common pathways involved in epigenetic 
dysfunction. Such a pathway could be related to the strong asso- 
ciations found between the H3K27me3 and H3K9me3 stem 
cell domains and aberrant DNA methylation. In general, there 
appear to be genomic locations which are highly susceptible to 
aberrant DNA methylation during toxicant-mediated malignant 
transformation. 

In the arsenical-transformed UROtsa models we observed 
only one potential region of agglomerative DMRs (Fig. S9), 
which occurs at the HOXB gene family cluster. The lack of 
agglomerative DNA methylation events in the arsenical-trans- 
formed UROtsa cells relative to CAsE-PE appears related to the 
epigenetic state of the parental UROtsa prior to arsenic exposure. 
The pathway of events that gave rise to the immortal UROtsa 
cell line induced pronounced epigenetic remodeling. This is 
apparent in the comparison of UROtsa with HUC, which shows 
many agglomerative DNA methylation events are present in the 
UROtsa cells prior to any arsenic exposure (Fig. S10). This find- 
ing is not unprecedented since others have also observed agglom- 
erative DNA methylation associated with catastrophic events 
such as overcoming stasis or immortalization. 14 This can help 
explain why arsenic did not induce as many agglomerative events 
in the arsenical-transformed UROtsa models compared with 
CAsE-PE, because genomic regions likely targeted for agglom- 
erative DNA methylation during arsenic exposure were already 
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Figure 5. Protocadherin G gene family cluster heatmap generated from the custom Agilent tiling MeDIP-on-Chip microarray data. The color of each 
tile in the heatmap represents the relative degree of methylation. Refseq genes, CpG islands and CAsE-PE (n = 2) hyper DMRs are shown below the 
heatmap in their corresponding locations. The small black bars mark the regions that were confirmed by the MassARRAY EpiTYPER assay. The smaller 
heatmaps were generated from the MassARRAY EpiTYPER data and the coloring of each tile represents the degree of methylation of each CpG unit 
from 0%to 100%. 
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Figure 6. Total number of hypermethylated probes from CAsE-PE 
(n = 2) that fall within H3K27me3 or H3K9me3 stem cell domains. 76% 
of the hypermethylated probes are within H3K27me3 domains (orange) 
when only 27% were expected. Eighteen percent of the hypermethyl- 
ated probes are within H3K9me3 domains (green) when only 4% are 
expected (***p < 1 x 10 ,6 , hypergeometric test). 



hypermethylated in the parental UROtsa. These observations 
also highlight RWPE-1 as a valuable model to study the epigen- 
etic mechanisms of toxicant-induced malignant transformation 
because its DNA methylation profile is relatively normal despite 
being an immortalized cell line. 



Several important similarities and connections exist between 
stem cells and cancer cells that support an epigenetic compari- 
son of these distinct cell types (reviewed in ref. 48). First, both 
exhibit self-renewal capability. Second, cancer cells often lose 
the terminally differentiated phenotype of their normal tissue 
counter-parts to resemble undifferentiated or stem cells. Finally, 
all normal differentiated cells and cancer cells are descendants 
of stem cells. These similarities and observations provide ratio- 
nale for a comparison of the distinct epigenetic states between 
normal embryonic stem cells and our malignantly transformed 
models and suggest that the epigenetic state of malignant cells 
could be related to or partially determined by that of their stem 
cell ancestors. 

As expected, we found that H3K27me3 stem cell domains 
significantly overlapped individual hyper DMRs that arise dur- 
ing toxicant-induced malignant transformation (Fig. 6). An 
association between H3K27me3 domains and aberrant DNA 
methylation in cancer has been described previously, which 
lends support to our data analysis because it conforms with the 
findings of several other studies. 40,42,43 When we compared the 
enrichment of H3K27me3 and H3K9me3 domains to CAsE-PE 
agglomerative DNA methylation only the H3K9me3 enriched 
regions had a striking overlap with the agglomerative DMRs 
(Fig. 4). This supports the link between H3K27me3 and aber- 
rant hypermethylation of individual genes, but in the case of 



www.landesbioscience.com 



Epigenetics 



1245 



agglomerative hypermethylation, H3K9me3 may be a more 
important player. 

The significant overlaps observed between stem cell H3K9me3 
enriched domains and the toxicant-induced agglomerative 
DMRs suggest that H3K9me3 modifications may be involved 
in directing aberrant DNA methylation to these regions. While 
the association of aberrant DNA methylation with H3K9me3 
in humans is rather new, normal mechanisms are known which 
link these two epigenetic modifications. For example, both are 
hallmarks of normal mammalian heterochromatin and the his- 
toric methyltransferases Suv39h and G9a have been shown to 
recruit DNA methyltransferases which can lead to spreading of 
DNA methylation in normal cells. 44 " 46 On the other hand, there 
must also be normal mechanisms leading to H3K9me3 that do 
not induce DNA methylation since many H3K9me3 domains in 
normal cells are not DNA methylated. 45 Perhaps during malig- 
nant transformation there are aberrant interactions between 
histone methyltransferases and DNA methyltransferases (e.g., 
Suv39h and DNMT3B) that lead to aberrant DNA methylation 
of H3K9me3 domains that are normally unmethylated. Once 
such a region contains aberrant DNA methylation, H3K9me3 
and its associated factors may induce further spreading of the 
aberrant methylation to result in an agglomerative region of 
aberrant DNA methylation. 45,46 

In conclusion, agglomerative DNA methylation is associated 
with toxicant-induced malignant transformation. Overall, this 
study suggests that there are common pathways of epigenetic dys- 
function in distinct models transformed by arsenic, cadmium or 
MNU, which target both individual genes as well as larger chro- 
mosomal regions that contain many genes. Individual DMRs 
are likely to occur in H3K27me3 domains while the agglomera- 
tive DNA methylation may be targeted to H3K9me3 domains. 
These associations provide further insights into the epigenetic 
mechanisms of toxicant-induced malignant transformation and 
carcinogenesis in general. 

Materials and Methods 

Cell culture and sample acquisition. UROtsa, URO-ASSC and 
URO-CDSC cells were contributed by Drs Donald and Mary 
Ann Sens. The UROtsa cell line was created from the urothe- 
lial cells of a 12-y-old female donor that were immortalized 
using a temperature sensitive SV40 large-T antigen construct. 25 
URO-ASSC cells were created through continuous exposure of 
UROtsa to 1 (xM arsenite for one year. 5 Similar to the arsenite- 
induced transformation, URO-CDSC cells were created by con- 
tinuous exposure of UROtsa to 1 |xM cadmium for one year. 5 
URO-MSC cell lines were contributed by Dr Jay Gandolfi and 
were created through continuous exposure of UROtsa cells to 
50 nM MMA m as previously described. 31 Each of the UROtsa 
cell lines was cultured as previously described. 34 The blackfoot 
transitional cell carcinoma lines (BFTC-905 and BFTC-909) 
were acquired from the German Collection of Microorganisms 
and Cell Cultures. 39 50 BFTC-905 and BFTC-909 were derived 
from urothelial carcinomas of individuals who lived in black- 
foot endemic regions of Taiwan where high level exposure to 



arsenic is common. BFTC-905 and BFTC-909 were cultured 
at 37°C with 5% C0 2 in DMEM with 10% FBS and 1% vol/ 
vol penicillin-streptomycin solution with media changes three 
times per week and sub-culturing before reaching confluence. 
The non-tumorigenic RWPE-1 cell line was derived from pros- 
tate epithelial cells that were infected with human papillomavi- 
rus 18 to induce immortalization. 26 RWPE-1 cells were cultured 
in K-SFM with 50 ixg/mL BPE, 5 ng/mL EGF, and 1% vol/ 
vol antibiotic-antimycotic solution. RWPE-1 cells were con- 
tinuously exposed to 5 |xM sodium arsenite for up to 33 weeks 
and are referred to as CAsE-PE cells after undergoing arsenite- 
induced malignant transformation. 6,35 RWPE-1 cells were also 
malignantly transformed by exposure to 10 |jlM cadmium for 
10 weeks and the resulting line is called cadmium-transformed 
prostate epithelial cells (CTPE). WPE1-NB26 is a malignant cell 
line that was created previously by exposure of RWPE-1 cells 
to 100 ixg/mL N-methyl-N-nitrosourea (MNU) for 4 cycles of 
1 h/passage, after which cells were selected by xenografting twice 
through nude mice and growth on soft agar. 51 WPE1-NB26 cells 
were cultured as previously described. 51 Normal human pros- 
tate epithelial cells (PrEC) cryopreserved at passage two were 
acquired from ATCC and cultured in PrEBM™ with added 
BulletKit™ (Lonza). Normal female human urothelial cells 
(HUC) were acquired from ScienCell™ Research Laboratories 
and cultured in Urothelial Cell Medium (ScienCell Research 
Laboratories). The prostate cancer cell lines PC3, LNCaP and 
DU145 were cultured in Iscove's Modified Dulbecco's Medium 
(IMDM), RMPI-1640 and IMDM respectively with 10% FBS 
and 1% vol/vol antibiotic-antimycotic solution at 37°C and 5% 
C0 2 . Short tandem repeat (STR) profiling was performed on 
each cell line either to validate its origin or confirm a unique 
DNA profile in the case where none has been published. DNA 
samples from six human urothelial carcinoma tumor biopsies 
were purchased from Origene Technologies. 

Nucleic acid isolation. Genomic DNA was isolated using 
the DNeasy Blood and Tissue Kit according to manufac- 
turer's protocol (Qiagen). The quantity of each sample was 
assessed using absorbance at 260 nm on the NanoDrop 1000 
Spectrophotometer. 

MeDIP coupled to microarray (MeDIP-on-Chip). MeDIPs 
were performed as previously described. 52 MeDIP and input 
samples were amplified, labeled, hybridized and washed as pre- 
viously described. 53 The following modifications were made for 
samples hybridized to the G4489A Human Promoter ChlP-on- 
Chip Microarray Set 244K designs 014706 and 014707 (Agilent 
Technologies). Two micrograms of each amplified MeDIP and 
input were labeled for each slide of the 2-slide set. After labeling, 
the MeDIPs and inputs were combined and then divided into 
two equal volumes before drying down to 187.5 u,L each. To each 
was added 12.5 (JtL human cot-1 DNA (1 mg/mL, Invitrogen 
Cat. No. 15279-011), 50 uX 10x Agilent blocking agent and 
250 ixL 2x Agilent hybridization buffer before incubation at 
95°C, 37°C and hybridization. Independent MeDIPs and inputs 
from identical or similar samples to those used on the Agilent 
G4489A Human Promoter ChlP-on-Chip Microarray Set 244K 
were amplified, labeled and hybridized to a custom Agilent tiling 
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microarray, design 024305, as previously described. 53 This micro- 
array covers 718 miRNA hairpin regions and includes extensive 
tiled coverage of all the HOX and protocadherin gene family 
clusters. 

MeDIP-on-Chip data analysis. Microarray data were 
imported to R 54 and normalized as previously described. 53 M val- 
ues (log2 ratios of immunoprecipitated and input channel) were 
used for further analysis as a measure of enrichment of a region 
centered on individual probes. The multidimensional scaling plot 
was generated using pairwise distances between samples that were 
derived from all M values and the function cmdscale in R. For 
calculation of DMRs, M values were analyzed in a sliding win- 
dow of 1,200 bp, with a step of one probe. DMRs were defined as 
regions of at least three consecutive probes less than 600 bp apart 
where the mean difference of ratios was at least 1.5 fold. A p value 
cut was selected to maintain the false discovery rate (FDR) < 1%. 
The FDR was determined by analysis of permutated data. 

Detection of agglomerates and regions enriched with stem 
cell domains. The Agilent G4489A promoter microarray cover- 
age was binned and annotated to individual genes. Each anno- 
tated bin was considered positive if it overlapped the domain of 
interest (DMR, H3K27me3 or H3K9me3 stem cell domain). 
Data were analyzed in windows of 49 consecutive bins with a step 
of one bin using the hypergeometric test for the difference of the 
expected and observed number of positive bins in each window. 
P values were adjusted for multiple testing. 

MassARRAY EpiTYPER™ DNA methylation analy- 
sis. EpiTYPER DNA methylation analysis was performed as 
previously described (Sequenom). 14,15 The EpiTYPER primer 
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